Is QGP really a liquid ? 
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Abstract 

Here we address the question regarding the nature of quark gluon plasma (QGP), whether it is a liquid 
or strongly coupled plasma (SCP), using two different phenomenological models, namely quasi-particle 
model (qQGP) and strongly coupled quark gluon plasma (SCQGP). First we compare these two models, 
both of which explains the results of lattice simulation of quantum chromodynamics, as a function of 
plasma parameter and conclude that the QGP is largely (T > 1.5T C ) SCQGP and only for T < 1.5T C it 
C* . ' ma y be a liquid. 
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quark gluon plasma. 



1 Introduction : 

The non- ideal behavior of QGP near T = T c , observed in lattice simulation (LGT) [T of quantum chromo- 
dynamics (QCD) and relativistic heavy ion collisions (RHICs) experiments, is still a highly debated issue in 
the study of QGP [2]. There are various QCD based approaches, perturbative [3], nonperturbative [4] and 
phcnomenological models to explain this problem. QCD based approaches, eventhough questionable due to 
strong coupling constant near T = T Cl was partially successful [3113]. At the same time, phenomenological 
models were successful to explain, but involves many fitting parameters. For example, largely studied qQGP 
models JSJUJ] has more than 3 parameters and recently introduced new qQGP 7, 8 J has still one parameter. 
Similarly, strongly coupled QGP [9lll0j which is inspired by similar studies in electrodynamic plasma, which 
is called SCQGP, has 2 fitting parameters. Both of these phenomenological models, qQGP and SCQGP, fit 
very well the LGT results, eventhough they are developed based on different philosophy. qQGP is based on 
the fact that the thermal properties of QGP may be understood by studying a system of non-interacting 
quasi-partons with effective mass related to collective properties of QGP, instead of a system of real partons 
with QCD interactions. Whereas in SCQGP, it is assumed that, just like in strongly coupled electrodynamic 
plasma (SCP), the basic interactions, QCD, is approximated by Coulomb interactions between partons and 
hence the equation of state (EoS) of SCP is modified for QGP. EoS of SCP was derived by partially analytic 
and partially numerical methods and parameterized as a function of the plasma parameter, T, which depends 
on the collective properties of the system. That is, T is defined as the ratio of average potential energy to 
the average kinetic energy. 

It should be noted that SCQGP [HI QH] is different from popularly known sQGP. sQGP means strongly 
interacting QGP in the sense that the coupling constant a s is not very small or weak and hence leads to 
non-perturbative effects like the existence of hadrons, mostly colored hadrons. It is claimed |Ilj that pressure 
due to these hadrons reproduces the pressure seen in LGT simulations. Of course, there is some confusion 
in using sQGP and SCQGP as discussed in [T5]. Later, Gelman, Shuryak and Zahed introduced another 
model cQGP (classical QGP) [13], based on molecular dynamics simulations, which is similar to SCP or 
SCQGP, but again claimed to be sQGP. Further, it is speculated that sQGP (or cQGP) might exhibit 
different phases like gas, liquid and solid as in the theory of SCP. Finally, above sQGP is now widely quoted 
as liquid QGP and based on the smallness of 2 ratio, it is also called perfect liquid, etc. Note that the EoS 



of cQGP, expanded in terms of T, to the lowest order, is very similar to one used in Ref. [HI HH] or SCP. 
Along with this work, Lin and Shuryak [14j also modelled this system by borrowing similar theory from 
super symmetry, string theory etc., regarding the evaluation of - so on and hot discussions on the bounds 
of -. But, it should be noted that gas, liquid and solid phases in a system with color charged partons is, in 
fact, SCQGP, analogous to SCP in plasma physics and not sQGP. Recently, Chernodub and Zakharov [T5] . 
and Liao and Shuryak [16] , independently, proposed a new picture to the existence of non- ideal nature of 
QGP near T = T c in terms of magnetic monopoles etc. and yet another idea to explain non-ideal behavior 
of QGP and it goes on. 

Here we address an interesting observation that the qQGP and SCQGP, seemingly two different models, 
both explains the LGT results very well. Of course, both the models are based on the collective properties 
of QGP, but the formulations are different. We compare the results of these two models and comment on 
the liquid nature of QGP. 

In section 2 we first briefly review SCQGP model, results and comment on the liquid nature of QGP. In 
section 3 a brief review of qQGP and results are presented. In section 4 we compare above two models as 
a function of the plasma parameter T. Results and conclusions on the comparison of two models and the 
predictions of the models on the nature of QGP is prsented in section 5. 

2 SCQGP: 

SCQGP was first proposed in 1999 [9] and it explained remarkably well the LGT results on gluon plasma 
and recently, it is also applied to flavored QGP and found to explain the LGT results very well [TU]. Here 
one modifies the EoS of SCP to SCQGP [TU] as, 

e(r) = — = i + -Uex(r) , (i) 

£sb 2.7 

where Esb is the Stefan-Boltzmann gas limit of QGP and u ex (T), is given by, 
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(3) 
i°° p = -0.898004 T + 0.96786 T 1/4 + 0.220703 r~ 1/4 - 0.86097 . (4) 



The excess energy u ex due to non-ideal interactions in above equation is obtained from SCP. SCP, a system of 
quasi-neutral charged system with Coulomb interaction, is extensively studied and the EoS is parameterized 
as a function of T and tested for a wide range of T, V < 180, and different systems like white dwarfs, dusty 
plasma etc. In QGP also we expect Coulomb like interactions, due to one-gluon-exchange, as in the case of 
hadron spectroscopy. Contrary to hadron spectroscopy, here in QGP, the confinement interaction between 
partons are neglected. Hence, SCQGP may also behave like SCP and T for QCD [9l [10] may be written as, 

<PE> _ 9cf^ 
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where r av may be taken as (jf^-) 1 ^ 3 with density n « 1.1 ajT 3 and aj = (16 + 21 nf/2)n 2 /90. We choose 
phenomenological running coupling constant as, 

(in 



(5) 



<**(T) = 



3(153 - 1971/) ln(2 ln(T/A T )) 
(33 - 2n/ ) m(T/A T ) V 1 ~ (33 - 2n f ) 2 ln(T/A T ) 



(G) 



Of course, many authors [TTl [TH1 US] included a factor 2 in the definition of T for QCD plasma to include the 
magnetic interactions, which amounts to a different value for the fitted parameter g c in our model. In our 
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Figure 1: Plot of A = (e — 3P)/T 4 as a function of T/T c for gluon plasma using the models SCQGP (A s 
qQGP (A,) and LGT results (symbols) [T]. 



model g c and At are the two parameters which may be different for different systems |10) like gluon plasma, 
2-flavor, 3-flavor and (2+l)-flavor QGP. Here we consider a simplest system, gluon plasma, and hence nj = 



in a/ and a s (T). To fit the LGT results [T], we used g c — 2.8 and to = At/T c = 0.5 and obtained a very 
good fit. As an example we plotted A = (e — 3P)/T 4 as shown in Fig. 1 (A s ) along with that of qQGP 
model and LGT results [T] and the resulting T(T) is plotted in Fig. 2. The plasma parameter T is about 0.5 
at T — 5T C and only close to T — T c , it rapidly grows to 1.8. If we take the results of SCP [20] that plasma 
is in liquid state for T > 1, QGP is not in liquid state except very close to T = T c , i.e., for T c < T < 1.5T C , in 
contrary to earlier claim in Ref. [TTJ [TH] [19] . Note that this conclusion follows from strongly coupled plasma 
formalism of QGP which fits LGT results. Even if we take a factor of 2 due to magnetic interactions our 
conclusions will not change, but the value of fitting parameter g c will be reduced by a factor of half. It is 
interesting to note that the plot given in Ref. [T7] for r(T), using qQGP, is similar to Fig. 2 if one scales it 
by a factor of half and T is less than 1 up to T — 1.5T C and increases rapidly to greater than 1 as T decreases. 
Probably the liquid-gas phase transition of QGP and confincment-deconfinement phase transition or cross 
over may be coinsiding and one needs a quantitative description of liquid-gas transition and the critical value 
of r (r c ) to fully understand the phenomena. Until then it is too drastic to speculate the existence of liquid 
and solid or crystal phase of QGP as in [T51[TB] and, from Fig. 2, the maximum value of T ss 1.8 at T = T c . 
It is interesting to note that recently Chernodub and Zakharov [T5] also have shown that QGP may be in 
liquid state only for T c < T < 2T C . 

Thus we see that nonrelativistic, classical strongly coupled plasma model with phenomenological mod- 
ifications to take into account of relativistic, flavor, color factors and quantum effects in terms of running 
coupling, fits well LGT data in the relevent temperature range, T c to 5 T c , and shows that QGP is strongly 
coupled plasma and may not be liquid. Of course, this model has a limitations that it is a classical model 
and hence for weak coupling limit it goes to Debye-Huckel results and not the perturbative QCD results. As 
shown in Ref. [5T], perturbative results also goes to Debye-Huckel term for classical limit. 

3 qQGP: 

As we discussed earlier, there are varieties of qQGP, starting from the work of Goloviznin and Satz 22J, 
Peshier et. al. [5] and many other authors [6] [3 [8] . Recently, we developed a new qQGP pj [8] which differs 
from earlier qQGP [5] [S] in the derivation of thermodynamic (TD) quantities from statistical mechanics 
(SM). We start from energy density a well defined in SM, and derive all other TD quantities using TD 
relations. Whereas in other models, one starts from pressure [5] and hence TD inconsistency and the need 




Figure 2: Plot of T as a function of T/T c for gluon plasma. 



of reformulation of SM [6] . In our qQGP model also there are different version [8] because of the way we 
choose the quasi-parton masses. There are varieties of qQGP models as well, depending on different forms 
of thermal masses and running coupling constants [6]. All models, by adjusting the free parameters, explain 
very well the LGT results. All these qQGP may be called classical. There is another class of qQGP, called 
dynamical qQGP, which follows from field theoretical approach. Again, one starts from pressure-partition 
function relation, but thermodynamically consistent because of the self-consistent $ derivable approximation 
[H [S3] . It reproduces leading order perturbative QCD results at weak coupling limit for gluon plasma and 
is successful in fitting the LGT results for T > 2.5T C [U[S3], but one needs to use phenomenological model 
with 2 or 3 parameters to fit the LGT results [24] for T < 2.hTr. 

Let us consider our simplest qQGP model where we take the thermal mass of gluons is equal to the 
plasma frequency and it fits the LGT results very well pj [8] . We start with the energy density pj [8] for 
gluon plasma, 

(7) 
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where e^ = \/k 2 + m% and /3 = ™. m g is the thermal mass which is taken to be equal to the plasma 



frequency 
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with the same a s (T) as given by Eq. ©. We assumed here that the whole thermal 



energy is used to excite quasi-particles and hence the vacuum energy is taken to be zero. As pointed out in 



Ref. [5], there is no TD inconsistency in our model with this assumption. From above equation, Eq. ([7]), we 

get, 
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Note that the pressure, in both the models, is obtained from the TD relation, 

dP 

e = T P 

dT ' 



(8) 



(9) 



on integration, with one integration constant which we fix to the LGT data at T — T c . The adjustable 
parameter is only Ay which we adjust such that we get the best fit to LGT results and A^/T c = .65. As an 
example we plotted A q in Fig. 1. This model is very successful to explain LGT results on all systems [7J|5]- 
Thus we see that this simple, single parameter qQGP model explains the LGT results as good as SCQGP in 
the relevent range of temperature T c to 5T C . But it has the limitation that it fails to reproduce exact leading 
order perturbative QCD results at weak coupling limit, evcnthough it leads to a term of the order of g 2 . 




Figure 3: Plot of e = — — as a function of T for gluon plasma. 



4 qQGP and SCQGP- Are they same ?: 



It is surprising why both models fit LGT results very well? Is it always possible to model SCQGP using 
qQGP and vise versa?. We see from Fig. 1 that they almost lie one over other. To address this question we 



plot e s and e q as a function of T in Fig. 3. This ratio, e s (r), for SCQGP is obtained from the Eq. |T]). The 
same quantity for qQGP, e q (T), may be obtained from Eq. l[5]). by expressing a s in plasma frequency in 
terms of T, using Eq. |5]). Note that in quasiparticle models the effects of interactions between real partons 
is replaced by the thermal mass of quasipartons. Therefore, the plasma frequency in our model includes 
the effects of interactions. Whereas in SCQGP the effect of interactions enter through plasma parameter T. 
Both plasma frequency and plasma parameter depends on a s , density and/or temperature and hence can 
be related. Both the models approximately match for T < 1, deviate rapidly for larger T. ft is interesting 
to note, from Fig. 2, that only very close to T — T c , T is greater than 1 and hence, except very near to 
T = T c , both models approximately match. Also note that the EoS of SCQGP, Fig. 3, shows a sign of phase 
transition, e s — > for some T, but qQGP does not show such features for any V and hence qQGP model 
may not model QGP liquid as claimed by [T7] . 

5 Results and Conclusions: 

We analyzed and compared two different phenomenological models of QGP, quasi-particle model and strongly 
coupled plasma model. The results are presented in Fig. f , where, as an example, A = (e — 3P)/T 4 of the 
models were compared with LGT results [Tj for gluon plasma and a remarkable good fit was obtained for 
both the models. Further, plots of the model almost lie one over other. In Fig. 3, the ratio e(F) = — — is 
plotted for both the model as a function of plasma parameter F. We found that SCQGP may be modelled by 
qQGP and vise versa for T < f , but models deviate rapidly as a function of T for larger T. By fitting the LGT 
results of gluon plasma, both the models predict that T is always less than 1, except for T c < T < 1.5T C , and 
hence QGP may be strongly coupled plasma state, rather than liquid or crystals state. For, T c < T < 1 .5T C , 
QGP may be in gas-liquid phase transition region which overlaps with deconfinement-confinement phase 
transition or cross over region which may contain many non-perturbative objects like colored and colorless 
hadrons |11] , monopoles [151 116) etc. which anyway make the region very blurred. Eventhough QGP is 
SCQGP one may use fluid theory to describe it's evolution as done in electrodynamic plasma where it is 
found that about 80 % properties may be explained. In this sense we can treat QGP as a liquid, as pointed 
out by Nagle [T^]. Similarly, the use of term sQGP may be a general one, including all models of strong 
interactions, QCD, but it may be useful to classify them with specific names like SCQGP [91 [10] or cQGP 
|13j . qQGP [551 13 HI HI [H] > colored and colorless hadron gas model [TT], supersymmetric model [II], monopole 



model [ini US] , so on and compare them. 
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